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Abstract In this paper we derive, by means of P-convergence, the shallow 
shell models starting from non linear three dimensional elasticity. We use the 
approach analogous to the one for shells and plates. We start from the min- 
imization formulation of the general three dimensional elastic body which is 
subjected to normal volume forces and free boundary conditions and do not 
presuppose any constitutional behavior. To derive the model we need to pro- 
pose how is the order of magnitudes of the external loads related to the thick- 
ness of the body h as well as the order of the "geometry" of the shallow shell. 
We analyze the situation when the external normal forces are of order /i", 
where a > 2. For a = 3 we obtain the Marguerre-von Karman model and for 
a > 3 the linearized Marguerre-von Karman model. For a € (2, 3) we are able 
to obtain only the lower bound for the i^-limit. This is analogous to the recent 
results for the ordinary shell models. 

Keywords Marguerre-von Karman model • Gamma convergence • shallow 
shell • asymptotic analysis 
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1 Introduction 

The study of thin structures is the subject of numerous works in the theory of 
elasticity. Many authors have proposed two-dimensional shell and plate models 
and we come to the problem of their justification. There is a vast literature on 
the subject of plates and shells and some also devoted to shallow shells (see [U 
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[2]). The expression shallow shell means that the curvature of the mean surface 
is also small with respect to the sizes of the mean surface. 

The derivation and justification of the lower dimensional models, equi- 
librium and dynamic, of rods, curved rods, plates, shells, shallow shells in 
linearized elasticity by using formal asymptotic expansion is well established 
(see [His] and the references therein). In all these approaches one starts from 
the equations of three-dimensional linearized elasticity and then via formal 
asymptotic expansion justify the lower dimensional models. One can also ob- 
tain the convergence results. Complete asymptotic (with higher order terms) 
for linear shallow shells and the influence of boundary conditions (boundary 
layer) is discussed in [4 . 

Formal asymptotic expansion is also applied to derive non linear models 
of plates and shells (see [Tl[3] and the references therein) , starting from three- 
dimensional isotropic elasticity (usually Saint- Venant-Kirchoff material). Hi- 
erarchy of the models is obtained, depending on the boundary conditions and 
the order of the external loads related to the thickness of the body h (see also 
[5] for plates). Formal asymptotic expansion is also applied to derive non linear 
shallow shell models (see [HEIIT]). It turns out that the asymptotic analysis of 
the equations of shallow shells closely follows that of a plate. If we assumed 
that the order of the normal external loads behaves like h^, we would ob- 
tain Marguerre-von Karman equations. Influence of different lateral boundary 
conditions on the model is also discussed (see [51[n])- 

However, formal asymptotic expansion does not provide us a convergence 
result. The first convergence result, in deriving lower dimensional models from 
three-dimensional non linear elasticity, is obtained applying r'- convergence, 
very powerful tool introduced by Degiorgi (see [TUIITT] ). Using /^-convergence, 
membrane plate and membrane shell models are obtained (see [121113] ). It is 
assumed that the external loads are of order h^. The obtained models are 
different from those ones obtained by the formal asymptotic expansion, at 
least for some specific deformations (compression). 

Recently, hierarchy of models of plates and shells is obtained via F- con- 
vergence (see P^[T^[TBlll7l[T5] l. Influence of the boundary conditions and the 
order and the type of the external loads is largely discussed (see [T51IT9] ). Let 
us mention that T- convergence results provide us the convergence of the global 
minimizers of the energy functional. Recently, compensated compactness ar- 
guments are used to obtain the convergence of the stationary points of the 
energy functional (see [301121] )■ 

Here we apply the tools developed for plates and ordinary shells to ob- 
tain shallow shell models by /^-convergence. It is assumed that we have free 
boundary conditions and that the external loads are normal forces with order 
a greater than 2. The main result is given in Theorem [Sj For the normal vol- 
ume forces of order a = 3 we obtain Marguerre-von Karman equations and 
for the order a > 3 we obtain linearized Marguerre-von Karman equations. 
Thus we have also justified linearized theory from three-dimensional non linear 
elasticity postulating only the smallness of the forces. We do not presuppose 
any constitutional behavior of the material. 
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In the situations when we have specific geometry of the shallow sheU (devel- 
opable surface) we are also able to obtain the model for a € (2,3). Otherwise, 
greater influence of the geometry of the shallow shell is expected (see Remark 

ED. 

Throughout the paper A or {A}^ denotes the closure of the set. By a 
domain we call a bounded open set with Lipschitz boundary. I denotes the 
identity matrix, by S0(3) we denote the rotations in M.^ and by so(3) the set 
of antisymmetric matrices 3 x 3. a;' stands for (xi, X2)- ei, 62, 63 are the vectors 
of the canonical base in R^. By id we denote the identity mapping id{x) — x. 
By V/i we denote V/i = Vei,e2 + ^^eg- — > denotes the strong convergence 
and ^ the weak convergence. We suppose that the Greek indices a, /3 take the 
values in the set {1,2} while the Latin indices i,j take the values in the set 
{1,2,3}. 

2 Setting up the problem 

We consider a three-dimensional elastic shell occupying in its reference con- 
figuration the set 0'^}', where ©''(I?''), fi''- ^ uj x {-h,h), w is a 
domain in and the mapping 0^ : {17''}^ R^^ is given by 

e'^ix'') = {xuX2,e''ixuX2)) +x'^a'^ixi,X2) 

for all x'^ — {x',x^) £ Sl^, where is a unit normal vector to the middle 
surface 0^{uj) of the sheU. We assume that 6'^ = where 9 € C^(a;), 

lim/^_s.o ~ 0, > 0. Using that assumption we can conclude that at each 
point of the surface ui the vector ag is given by 

"3 = («'r'^'(-/'5i0,-/''920,i), 

where 

a'' = {f''f\di9f + {f'')^\d29\^ + l. 

By inverse function theorem it can be easily seen that for h < ho small enough 
0'* is a diffcomorphism (the global injectivity can be proved by adapted 
compactness argument, see |3i Thm 3.1-1] for the ordinary shell). The following 
theorem is easy to prove and is a direct consequence of Theorem 3.3-1., page 
219, 220 in [T]. 

Theorem 1 Let the function 9^ he such that 

9'\xi,X2) = f''9{xi,X2), for all ixi,X2) € ui, 

where 9 G C^i^j) is independent of h. Then there exists h^ = hQ{9) > such 
that the Jacobian matrix's/ 0^ {x^) is invertihle for all x^ G and all h < Hq. 
Also there exists C > such that for h < ho we have 



detV0'' = l + {f'')^S''{x''), 



(2.1) 
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and 



V&\x'^) = I - /''C(x') + max{(/'')2, /i(/'')}R^ (x'^), 
(Ve^x"))-! = I + f^C{x') + nmx{(/'^)2, /i(/'')}R^(a;"), 

||(V6>'') - I||ioo(j7h.R3x3), ||(V6>'')^"' - I||i=o(r2'>;R3x3) < C/'*, 

1 V0''(a;',a;;j + /is)- V0''(a;',x?^) 



<C, 



where 



and 5^ : 



C(x') = 



5i6l(x') 

d2elx') 

-die{x') -d20ix') 

^3x3^ fc = 1, 2 are functions which satisfy 



(2.2) 
(2.3) 
(2.4) 

(2.5) 
(2.6) 



sup max \6'^{x'^)\ < Cq, sup max max |Rj^„(a;'')| < Cq, k = 1,2, 
o<h<hox'^en'^ o<h<ho x'^eo'^ 

for some constant Co > 0. 

Proof. It is easy to see (see also [Tj p. 220]) 



{V0\x''))n = 1 - ^x'^{a''{x'))-^^^2a\x')d,ieix') 

-dia'\x')di9ix')), 
{V&\x'')),2 = -^-^xl{a'\x')-^'\2<J\x')d^2e - d2a''{x')d,e{x')),{2.%) 



(Ve^a;'')) 



(V0''(a;'')) 



13 = ~f'\a\x'))-^'^d,e{x'), 
21 = -^xJ(a'^(x'))-3/'(2a"(x')ai2e(x') 
-d^a''{x')d2e{x')), 



{y0\x''))22 = 1 - ^x'^ia''{x')r'^^2a\x')d220{x') 

~d2a'\x')d2e{x')), 
iV0\x'^))23 - --f\a'^ix'))-'/^d29{x'), 

31 = f'^d.eix') ~ ^(a''(x'))-'/^d,a''{x'), 



{V&^{x^)) 



{V0'\x^)):,2 - f''d2e{x') ~ -^{a\x'))-^'^d2a\x'), 



{\/0\x^)h, ^ (a^{x')) 



(2.7) 
,(2.8) 
(2.9) 

(2.10) 



(2.11) 
(2.12) 

(2.13) 

(2.14) 
(2.15) 



Everything is already proved in '1 , except the relation (|2.5p . which is an easy 
consequence of the relations (|2.7p - (l2.15p . 
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The starting point of our analysis is the minimization problem for the 
shallow shell. The strain energy of the shallow shell is given by 

kHv) = f W{Vy{x))dx, 

where W : M^^'^ — [0, +00] is the stored energy density function. W is Borel 
measurable and, as in [MlfTSlfTBl . supposed to satisfy 

i) W is of class in a neighborhood of S0(3); 

ii) W is frame-indifferent, i.e., W^(F) — VF(RF) for every F e 'E?^'^ and 
Re S0(3); 

iii) W{Y) > CH'dist2(F,SO(3)), for some Cw > and all F e M.^""^, 
W{F) = if F e S0(3). 

By Q3 : R3><3 R we denote the quadratic form QsiF) ^ D'^W{T){F, F) and 
by Q2 ■ M^^^ M the quadratic form, 

Q2(G) = minQ3(G-|-a«)e3-|-e3«)a), (2.16) 

obtained by minimizing over the stretches in the X3 directions. Using ii) and 
iii) we conclude that both forms are positive semi-definite (and hence convex), 
equal to zero on antisymmetric matrices and depend only on the symmetric 
part of the variable matrix, i.e. we have 

Q3(G) =Q3(symG), g2(G) - Q2(sym G). (2.17) 

Also, from ii) and iii), we can conclude that both forms are positive definite 
(and hence strictly convex) on symmetric matrices. For the special case of 
isotropic elasticity we have 

Q,{F)=2^,\^^\' + X{trF)^ 

Q,{G) = 2H^^P + ^(tr G)^. (2.18) 

We suppose that the external loads are dead normal volume loads and thus 
we have that the total energy functional is given by 

j\y)=K'^iy)- [ fy,, 

where f:^ e L^{f2'').We suppose that the body is free at the boundary and the 
total energy functional is defined on the space W^''^{f2;M.^). Since the volume 
of physical domain decreases with the order 1 as ft. it is natural to look for 
the limit of the sequence j^J^- In fact, since the model crucially depends on 
the assumption how is the order of magnitudes of external loads related to the 
thickness of the body (see P^ITClfTB] '). we shall look for the limit of -^jiJ'^- 
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We shall analyze the situations when h "^E^ 0. For the applied forces we 
suppose 

^f:, - /3 in L\uj\ (2.19) 

where a > 2. The following questions have to be answered: For given order of 
external loads what is the order of the energy functional such that we have 
non trivial r'-limit? How does the limit functional look like? The main result is 
given by Theorem [5l where is answered to these questions for special cases of 
a. We take the special form of applied forces where the components of the force 
in ei and 62 direction vanish i.e. we suppose f \ = /2 = 0. For the analysis of 
the different situation see Remark [7] and Remark [5] Also we assume 

( fldx = 0. (2.20) 

This means that the force on the body is equal to and is a necessary condition, 
since it avoids the absence of a lower bound of the total energy functional 
arising from the trivial invariance y — ^ y + const. Firstly, we analyze how 
the order of the strain energy affects the limit displacement. Let us denote 

3 r'-convergence 

We shall need the following theorem which can be found in [T?]. 

Theorem 2 (on geometric rigidity) Let U C M™ be a bounded Lipschitz 
domain, m > 2. Then there exists a constant C{U) with the following property: 
for every v G W^''^{U;M."^) there is associated rotation R £ SO(m) such that 

II V-u - HWlhu) < C{U)\\ dist(Vt;, S0(m)|U2(c,). (3.1) 

The constant C{U) can be chosen uniformly for a family of domains which are 
Bilipschitz equivalent with controlled Lipschitz constants. The constant C{U) 
is invariant under dilatations. 

In the sequel we suppose ho > ^ (see Theorem [1]). If this was not the case, 
what follows could be easily adapted. Let us by P'^ : fi — > denote the 
map P''(x', X3) = {x' , 2/1x3). In the same way as in |T5i Theorem 10] (see also 
[IHl Lemma 8.1]) we can prove the following theorem. For the adaption we 
only need Theorem 13.11 and the facts that C{U) can be chosen uniformly for 
Bilipschitz equivalent domains and that the norms ||V6'''||, ||(V0'')~^|j are 
uniformly bounded on for h < ^. 

Theorem 3 Let a; C be a domain. Let be as above and let h < i. Let 
E = \ j dist2(Vy'\SO(3))dx. 
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Then there exist maps R'' : w ^ S0(3) and R'' : uj ^ R^'^^ , with |R''| < C, 

||(Vy'')o0''oF"-R"|U.(,,) <C\/:e, ||R''-R''|U.(^) <C\/^, (3.2) 

||V'R''iU2(„) < Ch-'VE, ||R'' - R''||loo(„) < Ch-'VE. (3.3) 
Moreover there exist a constant rotation Q'' e S0(3) such that 

||(Vy") o o - Q''|U2(,,) < CH-'Ve, (3.4) 

and 

IIR- Q''I|lp(.;) < C'p/i-i\/;B, Vp<oo. (3.5) 
Here all constants depend only on S (and onp where indicated) and f2 :— ]7^/^. 

Remark 1 By iWy^) o o P'' we have denoted Vy'* evaluated at the point 
0^{P''{x)). The construction of R'' and R'* is given in [IS]. Since S0(3) is 
a smooth manifold there exists a tubular neighborhood U of S0(3) such that 
the nearest point projection tt : U SO (3) is smooth. If E'^ < 6h^ then 
we always have R''(a;') G U. Hence the map R^ : S — )• S0(3), given by 
Rp(x') — 7r(R''(x')), is well defined and in Theorem [3] the map R can be 
replaced by the map R^. This is already noted in [15]. 

Remark 2 Since & is Bilipschitz map, it can easily be seen that the map y — > 
yo©'' is an isomorphism between the spaces W'^'^ f}''; R") and W'^'^{[2''; M™) 
(see e.g. [22]). 

To prove P-convergence result we need to prove the lower and the upper bound. 



3.1 Lower bound 

We need the following version of Korn's inequality which is proved in a stan- 
dard way by contradiction. 

Lemma 1 Let J7 C be a Lipschitz domain. Then there exists C{Q) > 
such that for an arbitrary u G W^^'^(I2;K^) we have 

l|M||wi.2(r2;R2) < C(X2)(|| symVM||L2(^2;R2) + I / udx\ + \ / {d2Ui- diU2)dx\). 

(3.6) 

Lemma 2 Let y'' e VF^'^ (!?''; R^) be such that 

\[ dist2(Vy'',SO(3))dx < Ci;'', (3.7) 
lim h-^E'' = 0. (3.8) 
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Let US also take = max{/i, /i E'^}. Then there exists maps R'' G 
W'^''^{uj, SO (3)) and constants R'^ G S0(3), c'' e such that 

t:^ (R'')V-c"- 
and the in-plane and the out-of-plane displacements 



U'\x') := 
V'^ix') := 



(.,3:3) - x')dx'i, 



'^'((y3^o0'^oP")(.,x-3)-m))rfx3 

-1/2 



satisfy 



L2(r2) 



(3.9) 



(3.10) 



IIR'' - IIIlpM < Cph-^VEh yp < 00, ||V'R''|U2(„) < Ch-'VEh. (3.11) 

Moreover every subsequence (not relabeled) has its subsequence (also not rela- 
beled) such that 



fEh 



V in W^'^iu:), v G W^''^{oj), 



u 



. fh"" 1 
mm I 



u m W^1'2(w;M2), 



fEh 



(R'*-I)^A m i«(cj;R^^'^), Vg < 00, 



((Vy'') o 0'' o P'' - I) ^ A inL^{n;l 

d3A = 0, A e PF1'2(w;R3x3)^ 

A = 63 (Ki V'w - V'u (g) 63, 
-^sym(R" - I) ^ ^ L'in;R^-') 



)3x3n 



(3.12) 

(3.13) 
(3.14) 

(3.15) 

(3.16) 
(3.17) 

(3.18) 



Proof. We shall follow the proof of Lemma 13 in [TS]. Estimates (13.101) and 
p. lip follow immediately from Theorem [3] and Remark [T] since one can choose 
R'' so that (1231) holds with Q = I. Using (jXTUl) and ((XTT|) we conclude that 

II (Vy") o o - < Ch-^y/E^. (3.19) 

For adapting the proof to the proof of Lemma 13 in [TS] it is essential to see 

(Vj/'*) oO^oP'' ^ (V(y'' o o P'')((V0'')-i o P'') 

= Vhiy'" o&'^o P''){{V0^')-^ o P''). (3.20) 
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From (IX^ it follows 

((Vy'') o e'' o P'')((V6>'') o P'^) = Vh{y" o o F''). (3.21) 
Using the fact that \\V0^\\i^a^(^^h) is bounded, p.l9p . p.21[) we conclude that 

||V,.(y'' o&'^oP'^- 0'^ o P'')|U2(^) < Ch-^VE^ (3.22) 

By applying additional constant in-plane rotation of order h^^\/E'^ to y'* and 
R'' we may assume in addition to p. 101) and (13. lip that 

[ {{d2y\) o 6>'' o P'' - {divl) ° 0'' ° = 0. (3.23) 

J n 

By choosing c'' suitably we may also assume that 

[ (y'' o o P'' - 6>'' o P'')dx = (3.24) 
J n 

From (I3.20p it can be easily seen that the following identity is valid 

{Vy^) o 0'^ o P'* - I = (V,,(y'' o 6>'' o p'' - ©'^ o P''))((V0'')-i o p'' - I) 
+(Vft(y'' o o P'' - o P'*)). (3.25) 

Using ([23]), dSmi), (12231), ((3:25)) we conclude that 

I I {d2iy^o0''oP''-0'loP'')-diiy2o&''oP''-0'^oP''))dx\ < Cfh-^^E^K 
Jn 

(3.26) 

Let us define A'' = (h/V E'^)(R^ — I). From (j3.1ip we get for a subsequence 
A''^A in M^1'2(^.r3x3) 

Using the Sobolev embedding we deduce p.l4p . Using p. 101) we deduce p.lSp . 
Since (R'')'^R'' = I we have A'' -I- (A'')^ = -(V^//i)(A'*)'^A'*. Hence A + 
A^ = and after multiplication by ^^^= we obtain p.lSp from the strong 

convergence of A''. From (12. 4p . p.22p and p.25p we conclude that 

-^V?,(y'o0''oP''-0'' oP'M^ A in L2(f2:M3x3)^ (3 27) 

Ve^ 

From (I3.24p . p.27p and the Poincare inequality we conclude the convergence 
in p. 121) . Moreover we have diV = A^i for i = 1,2. Hence v S VF^'^ since 
A e W^'^. Since A is skew-symmetric we immediately have A31 = div, A32 — 
d2V. If there exists constant C > such that h"*" < CE'^Vh < ho, we could 
multiply p.lOp with to conclude for some C > 

|||^ sym((Vy^) o ©'^ o P" - I) - |^ sym(R'' - I)|U.(^) < C. (3.28) 
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Using (j3.18l) we obtain that there exists C > such that 



1-^ sym((Vy'^) o 0- o - < C. 



(3.29) 



Using the identity (|3.25p we conclude that there exists C > such that 



sym(V,(y'^ o o - o P'^)|U.(,,) < C. 



(3.30) 



From that we conclude that \\^ sym V'M''||i2(^') is bounded. Using Lenima[TJ 
p.24p and p.26p we have the convergence p.l3p . In the situation where there 
exists subsequence such that liin;i_j.o h^^E'^ = we divide p.lOp with to 
conclude 



j= sym((Vy") o ©'^ o - I) - ^ syni(R'^ - l)\\mn) < C. (3.31) 



1 h'^ 

Since for that sequence we have < C-^ we conclude that 

II syni((Vy'') o O'^ o P^ - < C. 

In the same way as before we conclude that |1 SYmVu^\\i^2(^^-^ is bounded 
from which, again using Lemma [U p.24p and (I3.26p . we conclude convergence 
p.l3p . It remains to conclude A12 = 0. But this is easy, from (j3.13p and (j3.27p . 
using the fact that min{ , -^ttI-^^x" ~^ ^ as /i — 0. 

Lemma 3 Under the assumptions of Lemma\^we can also conclude that 



in (-^, -L=) [(y„ o 0'^ o P'^ - o P'^) + VE^x^d^v 



mm I 
h 



/Eh, 

(yg OO^OP^^- &l O P") 



all m W^'^{n). 

Proof. By the Korn's inequality and the relations p.27p . p.30p it is enough 
to prove that 



1 



o„ :— mm 



VP'*' y/Eh 



(y„ o 0" o P" - o P'* - [/^) + ^x^do.v 



fEh 



iVs o 0'' o P'' - 0^ opf"- V'') 0, 



^0, 
(3.32) 
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both in L^{fl). We shall just consider the situation h^^E'^ — )■ oo or h^^E'^ 
1. Using Poincare the inequality on each segment {a;'} x [— i, i], the Cauchy 
Schwartz inequality and p.27p we conclude 

J UJ 



< C- 



fEh 



£^''A„3||i2({^,}x[-i,i])C^2;' 



-A. 



a3 



L2(J7) 

(3.33) 



\o3\\L^n) 



<C^f dsiy^, o 0'^ o P^){x\ 

-a3(0^F")(x',.) 



< Ch- 



ids{y^,o0'^oP'^)-^d3{0'^oP'^) 



dx' 



L2(f2) 



0. (3.34) 



Lemma 4 Consider : ^ : w ^ S0(3) and E^ > Q and set 

1 



u := mm 



h r^l^ 



'^'\fy'lo0''oP'^\. , 



{{y'io0'^oP'^){.,x,)~f'^9{.))dx3. 



VE'"- J -1/2 

Suppose that we have a subsequence of y^ such that 

lim h-'^E'' ^ 0, 

/i-i-O 



and 



G" := 



(R'')'^((Vj/^') o 0^' o P'') I 

7w- 



(3.35) 
(3.36) 



\\{Vy^) o 0" o - n^Wmn) < C^, 
u^^u m 1^^'2(w;R2), v'' v inW^^^iuj), veW^'^iuj). (3.37) 
Then 

h 



(R'' - I) ^ A = 63 (g) V'u - V'z; (g) 63 m L^{u};R^''^), (3.38) 



G m i^(l2;R-*^-^). (3.39) 



and the 2x2 sub-matrix G" given by G^^ = G^^ /or 1 < ct, < 2 satisfies 
G"{x',x3) = Go{x')+x3Gi{x'), (3.40) 
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where 

Gi = -{V'fv. (3.41) 

Moreover 

V'u + (V'tt)'^ + V'w (g) V'l; + Vv ® Ve + Ve ® Vv = O, if h-^E'' oo, 

(3.42) 

symGo = i(V'M+(V''u)^+V'u0V'i;+V'u0V'6l+V'6'0V''y), ifR-^E^ 1, 

(3.43) 

symGo = ^(V'tt + (V'tt)'^ + V'w«)V'e' + V'6l®V'u), if h'^E^ ^ 0. (3.44) 

Proof. We follow the proof of Lemma 15 in [15]. We first assume (I3.38P 
(for the special sequence coming from Lemma [2] we know this anyhow) and 
establish the main assertion, namely the representation formula for G. Using 
the identity 2 sym(Q - I) = -(Q - I)^(Q - I) which holds for ah Q e S0(3) 
we immediately deduce from (I3.38P that 



Eh 



sym(R'^ - I) 



^(y'v^Vv + \V'vfe3(x)e3^ in L\uj). (3.45) 



By the assumption G'* is bounded in L^, thus a subsequence converges weakly. 
To show that the limit matrix G" is affine in we consider the difference 
quotients 

Il''{x',X3) = .s-'[G''{x',X3 + s)- G'^ix^x^)]. (3.46) 

By multiplying the definition of G'' with R'' and using (|3.2ip we obtain for 
{1,2} 

(R''H'')„^ = ^=[V„(?/'' o 0'^ o P>-){x', x3 + s)- 

-yh[y''oe''oph){x',Xi)U 

^ [((Vy")o0"oF")(a;',:E3 + s)- 
\-{{^&'')oP''){x',x^ + s) 



fEh 



a/3 



((Vy'^)o0"oP")(a;',X3)- 
l-{{V&'')oph){x',x^) 
{Vh{y'' o 0'' oph-0'^o P'')ix',X3 + s) - 
-V„(y'' o 0'^ o P'' - 0'^ o P''){x',X3) 



fEh 



((Vy'')o0''oP")(a;',X3 + ,s) 
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-R^x')) (l - ((V0'') o P^){x', X3 + s)) 
((Vy")o0''oP")(a;',X3) 
-R''(:r'))(l-((V0')oP'')(x',X3) 
(R''(a:') - I)((V0'') o P''){x',X3 + s 
^iiV0'')oP>^)ix',x,) 



Q/3 



(3.47) 



For the first term in p.47p we obtain 



^=[(V„(?/" o 0'' o P'^ - 0'' o P''){x', xs + s) 

-y^iy'^ o 0" op'^-e'^o p^){x\ xMcp 
= -^df, (- r (y" ° 0' ° P" 0' o 

VS'' Jo h 
From p.36p and (|3.38p we conclude 
h 



fEh 



{{yy^) o 0'' o P'^ - I) -> A in l2(I?;R='^3). 



(3.48) 



(3.49) 



From that, using that ||V0'' — IH^oofj-jh) — > and p.2ip . we easily conclude 
that 



Wh{y^ o0^o P'' - 0'=- oP^)^ A in L2(I2;M3x3)^ 



(3.50) 



Using p.50p we conclude that the first term in p.47p converges weakly in 
W-^'^{oj X (-1, 1 - s)) to A„3,;3(a;') = -w a/3(x'). By using (EH), (I23|), (1^:^ 
and (j3.38[) we conclude that the all the rest converges to in P^(ajx (— 1, 1 — s)). 
Since we have R'' — ^ I boundedly a.e. and H'* ^ H in we thus obtain 
Ha/aix' jXs) = —v^ap(x'). From that we conclude that G" is affine in 0:3 and 
that Gi has the form given in lemma. In order to prove formula for Gq it 
suffices to study 



G„^(x') 



G''(x',X3)(ix3. 



We have for a,/3 e {1,2} 



{G'')^p{x',x:,) 



((V'y") o 0" o ph - I)„^ 



fEh 



+ 



VPh 



, _ ^T{Vj/^)o0h_oPh_^ 

VPh 



(3.51) 
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In the case h~^E'^ 1 using ((^ , the identity ((X^ . (1530)) and the con- 
vergence of we have 



sym| 



1/2 ((V'y'^) o o P'^ - I) 



-1/2 



c^a^sla/J ^ sym[V'M + AC]a/3, (3.52) 



in L'^iuj). Using the convergence p.45p we easily conclude 

r A^i 
(symGo)Q/}^ sym V'm + syni(AC) - — in ^^(a;). (3.53) 

L 2 J a/3 

From p.53p we easily deduce p.43p . To derive p.42p let us observe that in 
the case h-'^E'' oo we have = h-^^/W. Multiplying ((33T|) by h?/^fW- 
(— ?► 0) we conclude 

-^[sym j^'\{Vv'^)o0'^oP'^-\)dx^]^p-sym[i^{K'^~l)U^OmL\io). 

(3.54) 

Using again (|2.2p , the identity p.25p and the convergence of u'* we have 



7,2 /■1/2 

-^[sym / ((V'y") o o P"- - l)dx:iU ^ sym[V'« + AC]„^ in L\uj). 

^ J-l/2 

(3.55) 

Using the convergence (I3.45P we conclude p. 421) . p. 441) can be concluded in the 
similar way from (13.511) . It remains to prove p.38p . Since R'' is independent 
of X3, we have for i, j e {1, 2, 3} 



(I 



i-l/2 J-l/2 



1/2 

-1/2 



(3.56) 

By using (j2.4L the relation p.25p and the convergence p.37p we conclude for 
a,/3e{l,2} 



II / ((Vy")o0'^oP"-I),^da:3||i.(,) < C{h-^^/W^f' + lTla.^{^\^}). 

J~l/2 " 

(3.57) 

From p.56p and p. 571) . by using p. 361) . we conclude 

h 



fEh 



(R'' - I)a^ ^ in L^{io) 



(3.58) 



In the same way, by using (|2.4p , the relation p.25p , p.36p and the convergence 
(|537|) . we conclude for ^ G {1, 2} 

II {{^y^) ° 0" ° - I)3/3rfX3 - (/^-^V^)a/3Z'|U2M 

J-l/2 

< Cih-^^/E^f^ + max{/E^, :^}). 



Shallow shell models by /^-convergence 



15 



In the same way as for p.58p we conclude 
h 



/Eh 



(3.59) 



Using the fact that R'' takes the values in S0(3) we deduce that ||Rj^3||i2((^-) < 
CV E'^/h. To get control on R33 we use the fact that for Q € S0(3) we have 

2 

|l-Q33| = |detQ-Q33|<C J2 l(Q-I)a^l+C(|Ql3Q3l| + |Q23Q32|). 

a, 13=1 

From this, using the generalized convergence theorem (with convergent ma- 
jorant rather than constant majorant), we easily deduce that {h/VW^CR^^^ — 
1) — )■ in L'^{uj). To control Rj'3 we use the fact that the first and third row 
of R'' are orthogonal. This yields 



|Rj3 + R§i|<C(|Rji-l| 



1 1^33 



|R'l2l), 



and together with p. 581) . p.59p and the convergence of R33 this gives the 
desired convergence for Ri^. The same argument applies to R23 and this 
finishes the proof. 

By /^^"-^ and /^^^^^ we denote the functionals 



/^'^"^(m,v)= f (lQ2{l{Vu+{V'uf + Vv(g>\/'v + Vv(g,\/'0 

(3.60) 
Vv) 



(3.61) 



/^^''^(m,v) = J (lQ2(i(V'M+(V'M)^ + V'v®V'0 + V'6 

+ ^Q2iiy'rv))dx', 

defined on the space P^-^'^(w;R^) x W^'^{uj). 

Remark 3 The term V'M+(V'M)^+V'w(g)V'w+V'-y(8)V'6'+V'6'(8)V' V measures 
the change of the metric tensor (which is of the second order), while the term 
(V')^w measures the curvature tensor change (which is of the first order) of 

' x' + {f^fu{x'y 

f^{v + e){x') 



the deformation 1^9(2;') 



with respect to the deformation 



ipo{x') 



f'Oix') 



(vV)^vV 



. To see this let us calculate 



(/'^)2V'm 

f''V'{v + 9){x') 



(/'')2V'm 

f''V'{v + 9){x') 



= I + {f'')\Vu + {Vuf + V{v + 0)(g> V'(« + 0)) + Oii/')^) 
= (VVo)^VVo + {f'^fiy'u + {Vuf + Vv (g) Vv + V'v (g) V'6 

+V'9(gVv) + 0{{f'')^). 
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The same explanation goes in the situations when h~*E^ — when the term 
Vv (S) V'ti disappears. To analyze the change of the curvature we need to 
calculate the normals 

- 63 - {f")di{v + e)e, - {f'')d2iv + e)e2 + 0((/'')2), 
n^o - 63 - {f)diee^ - {f")d29e2 + 0{{f''f). 

Now we have for arbitrary a, (3 £ {1,2} 

da^V ■ = f\d^p{v + 9)) + 0{{ff) 

Corollary 1 Let ,G^\G,G" .G^.Gi he as in Lemma\4\ 

Then we have the following semi- continuity results. 

i) //lim/j_j.o h^'^E'^ — oo then 

liminf-|^/'^(y'^)> f l^Q,{[Vrv)dx' . (3.62) 
a) Iflim.h^oh-^E^' = 1 then 

liminf ^/'Vy'') > /*^^-^(m,w). 
Hi) IfMuih^o h^'^E^ — then 

liminf 4r^''(y'') > I^^^"^ (u,v). 

Proof. We shall use the truncation, the Taylor expansion and the weak semi- 
continuity argument as in the proof of Corollary 16 in [15]. Let m : [0, oo) — >■ 
[0,oo) denote a modulus of continuity of D'^W near the identity and consider 
the good set Qh '■= {x € : |G^(a;)| < h~^}. Its characteristic function Xh is 
bounded and satisfies Xh — > 1 a.e. in fl. Thus we have XhG'^ G in L'^{f2). 
By Taylor expansion 

■^XhWil + /e^G'^) > iQsiXhG'^) - m(h-'^)\G'^f. (3.63) 

Using p. 351) and the boundedness of sequence G'' in L^{f2; M.^^^) we conclude 

liminf 4r-^''(y'') 

= liminf^ / W{{I{!'f{{Vy^)oe^'oP^))dx 
/i-i-o E Jq 

f Q3{XhG^)dx + ^ f {l-Xh)WiiVy'^)o0'^oP>^))dx 
i / Q^{G)dx > i / Q2{G")dx. (3.64) 

J ^2 J f2 



> 



> 

- 2 
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Here we have used the fact that Qa is a positive semi-definite quadratic form 
and therefore the functional v f^^ Qz(v) is weakly lower semi-continuous in 
1/2. Now by ((OO)) we have 

.1/2 ^ 

/ Q2(G")(x',X3)dx3 =Q2(Go(a:')) + -rQ2(Gi(:E'))- (3-65) 
Together with p.4ip . p.43p . (|3.44l) this implies the claim of the corollary. 



3.2 Upper bound 

Theorem 4 (optimality of lower bound) // h~^E^ — > 1 and if v € 

W'^''^{ijj), u e W^-'^iu:;^?) then there exists a sequence G W^''^{Q^]M.^) 
such that 

(Vy'') o 0'' o P'' ^ I m L^[2; M^xS)^ (3 66) 

and for , defined by kS. (where y should be replaced by y) convergence 
\3.12]) - fi3.1SI\) are valid and 



\^m^l\y^)^l''^^{u,v) 



(3.67) 



Ifh-'^E^ andifv ^ W'^-^{<jj), u £ W^-'^{oj;M.'^) then there exists ij^ such 
that the convergence IKWjl . iTO^j - iTO^) hold and 



1 



lim -=^r(y 
h^o E^ ^ " 



tLMvK 



iu,v). 



(3.68) 



Proof. Let us first analyze the situation when h "^E'^ — > 1 and assume that 
u, V are smooth. Then we define 



y'\0''ix',x'^)) = 0-ix',X3 
E'' 



^uix') 

^v{x'\ 

1 



/Eh 



'divix'Y 

4 I d2v{x') I + 



+ 7I^-3'^°(^')+2 h 



{xlfd\x'), 



(3.69) 



where do, di : a; M'' are going to be chosen later. The convergence (13. 66^ as 
well the convergences p.l2p - p.l3p can easily seen to be valid for this sequence. 
We also have 








V'v 






V'm 












/Eh 



x^d^ ® 63 



(3.70) 
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From (I3.70p by using (j2.3p we conclude 



h 










V'v 







3' 



fEh 



63 



(3.71) 



Using the identities (I + A)^(I + A) =1 + 2 sym A + A^A and (63 ® a' - 
a' ® 63)^(63 ®a'~a'® 63) ^ a' ® a' + |a'pe3 ® 63 for a' e we obtain 

(Vy'M^fVy'') = 1 + -^[2sym(V'tt + V'w® V'6i) + V'w® V't; + 

^/e^ E^ 

+ |V'wpe3 (g) 63] - 2—-x'i,{Vfv + 2— sym(d° ® 63) + 
ft /i-^ 



fEh 



x'i, sym(rf^ (g) 63) + 0{h^). 



(3.72) 



Taking the square root and using the frame indifference of W and the Taylor 
expansion we get 



(S'^)- V(Vy'') = {E^)-^W{[{Wy^fWy^f'^) 
where 



-Q3(A + a;3B), (3.73) 



A = sym(V'M + V'w ® + ^Vi; ® V'u + i|V'wpe3 g) 63 + sym(d° ® 63), 
B = -(V')^w + sym(d^(8)e3). 

For a symmetric 2x2 matrix A" let c = CA" £ M.^ denote the (unique) vector 
which realizes the minimum in the definition of Q2. i.e. 

Q2(A") = Q3(A" + c ® 63 + 63 «) c). 

Since Q3 is positive definite on symmetric matrices, c is uniquely determined 
and the map C is linear. We now take 

1 



^|V't;pe3 + C{sym{V'u + V'v (g) V'O) + V'v (g) V'v) 



= -2C{{V'fv). 

This finishes the proof of theorem in the situation h~^E^ — 1 and smooth tt, v. 
For general m, v it suffices to consider smooth approximations m'', u'', d"''', d}'^ 
of M,v in W^i'2(^) i Vl^2,2^^.]j2-)_ choose it'' e C°°(w; M^)^ ^ 

C°°(w) such that ||it''-M||vi/i,2 < ft and ||t;''-w||n.2.2 < ft and {u^ ,v^)- 



tMvK I 



{u, v)\ < ft. For u , V we choose y ' such that 



TMvK/„.h „,h 



(^t^^;'*)| < ft. 
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Then we have the claim. 

In the situation when h^^E^ — 0, for smooth tt,w, we define 

( ^fWi.iT' A a/p* ( divix'Y 

+^x^,dP{x') + l^{x1fd\x'). (3.74) 

In the same way as in the situation h~'^E^ — )• 1 we can define d} such that 
(I3.68P is satisfied (the term V't;(8) V'w disappears because it is of the order 
for which — is satisfied. For non-smooth u^v the argument is the 

same as before. 

Remark 4 The equality (I3.42p can be written in the form 

Vu + (Vuf = ve ® ve ~ V{v + e)® V{v + e). (3.75) 

The left hand side is the symmetrized gradient of a W^'^ function. It is known 
fact that, if lo is simply connected, a map e : w — >■ Kg^^ is the symmetrized 
gradient of a Vl^^'^(w;]R^) function it, i.e. 



if and only if 



2e = (V'm)^ + V'tt, (3.76) 



922611 + ^11622 - 2^12612 = 0, (3.77) 



in the sense of distributions. On the other hand it is easily seen that this 
condition on the right hand side implies (see Proposition 30 in 16 ) 

det(V')^(w + 6*) = det(V')^6'. (3.78) 

For smooth v this means that the graphs of the functions x' ^ {v + 6){x'), 
x' — > 0(x'^ have equal Gauss-Kronecker curvature at each point which is a 
necessary condition for the existence of an exact isometry between these two 
surfaces. 

Remark 5 Situation h~^E^ — )■ and lim/,,_j.o h~'^E^ -> oo remains uncovered 
for the upper bound. There are several reasons for that and they are similar to 
the one observed in [T^[T51[^[^ . In the situation h~^E^ — >• we need some 
additional regularity results (e.g. v G W^'°°{uj)). This can be concluded when 
the graph of 9 is developable surface i.e. when we have det(V')^^ — (see 
[15]). The situation h~^E^ > is more complicated. For given v G W'^''^{uj) 
we would like to construct the (exact or higher precision-see |24j ) isometry 
(from the graph of the function x' — f^9{x') to the graph of the function 
x' ^ f^{y + e){x')) of the form 



eh\2 
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The condition (I3.78P is a necessary condition for the existence of isometry, but 
not sufficient unless we have special situation det^(V'0) = K. In the situation 
K — one can construct an exact isometry in the similar way as in [15] (see 
also [17]). This can be done under additional (mild) hypothesis that there 
exists e > such that h~^'^^'^'> E'^ 0. The construction of an isometry would 
go by these steps: 

1. it can be seen that there exists the (^^(w) isometry between ui and the 
graph of f^6{u]) of the form 

where (j)h G C^(w) (see the proof of Theorem 25 in |TS])- 

2. Using this isometry and the one between w and the graph of f'^{v + 0){uj), 
one can easily construct an isometry between the graph of f'''9{uj) and the 
graph of f'^{v + 9){uj). 

In the general case, when 9{lo) is not developable surface, in the situations 
h~'^E^ — )• and h~^E^ — oo, stronger influence of the geometry of 9{uj) on 
the model is expected (see [H]). 



3.3 Convergence theorem 

Let fl e L^{Q^-M.) be given with the property 

/ f1 = Q, -^f.oe'^oP'^ ^ f,mL\n;R). (3.80) 

Jr}'! hVE"- 

Let m'' be the maximized action of force over all rotations of 1?'', 

m'' = max / f'i{x){Qx)3dx, (3.81) 

QGSO(3) 

and define 

M^iQe S0(3); r(Q) < +c5o}, (3.82) 

to be the effective domain of the following relaxation functional r : S0(3) — 

[0,+oo] 

r(Q) = minjliminf ^(m'^-y" /J(a;)(Q''x)3dx) ; Q'' G S0(3),Q'' ^ q}. 

" (3.83) 
The set A4 identifies the candidates for large rotation the body would perform 
to reduce its energy (see [IH1I23])- 
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Remark 6 The set M is introduced in |23) . We can conclude (using the uniform 
convergence of the functions Q — ^ 1^2 In^ f3{x)iQx)3dx) that ^^^^^ ~^ ^ 
where 

m= max f (f f3{x'.,X3)dx-:^{Ql(^ \)3dx' 

= max [( [ f:i{x' ,X3)dx3]{qiXi + q2X2)dx' . 

From that and the definition of M. and r it can be concluded that C A^o 
where 

= {Q € S0(3) : y ( y'J^ S3{x\x3)dx^ (Q {^'^)3dx' = m}. 

Since r is lower semi-continuous it can be concluded that is a nonempty 
closed subset of A^o- Since we also have the change of geometry as /i — >■ it can 
not be concluded that Al = A^n even under the condition — }= -f^oQ^ oP^ = 

/g for every h > 0. Let us mention that for every Q € A^o 3 A1 we have the 
following equality 

/ (/'^^//3(a;',a;3)rfa;3)(QF(^^'^)3dx' = 0, VF e so(3) (3.84) 

The equality p.84p is the consequence of the fact that the differential vanishes 
at the extreme points and the fact that so(3) is tangential to S0(3). The 
equality p. 841) is the balance of momentum. Since the "dead loads" are 
given on the reference configuration, the shell adjusts its deformation to satisfy 
the balance of momentum. 

To the total energy functional (defined on the space 14^^'^(i7''; M"^)) we add the 
constant and redefine 

J'iy") = I'-iy") 1 flyt + ^- (3-85) 

The following theorem is the main result. 

Theorem 5 (/^-convergence) Assume h^^E^ — > 1 or h^*E^ 0. Let us 
suppose that /J G L^(/2'';R) is given and satisfies Ii3.80\) . Then: 

1. There exists C > such that for every h > we have 

> inf |-i^j''(y'^);y'' e W^'^{n'';R^)^ > ~C. (3.86) 
^- //y'' G W^''^[fl^;^^ is a minimizing sequence of -^J^, that is 

then we have that there exists R,'* G SO (3), c'' G R such that the sequence 
(R,'', y'') has its subsequence (also not relabeled) with the following property: 
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i) For the sequence := (R'')^y'* - c'* and U'\ defined by iTO)] the 
following is valid 



/ 1 \ 

m'' := min (— ^= ^ it weafc/?/ m M^i'2(w;R^) 



^/so, any accumulation point R o/ i/ie sequence R'' belongs to A4. More- 
over if h^'^E^ — !■ 1 t/ien any accumulation point (tt,w,R.) of the sequence 
(m'', w'', R'') minimizes the functional 



J"(m,«,R) = R33 y (^j '^^J^dx^)v{x')dx' + r{K), (3 



w/iere J*^"-^ is rfe/inerf in ifX^Z!)] anrf r is defined in ^MW- U h-^E^ ^ 0_ 
then we have that any accumulation point (u, v, R) of the sequence («'', v'^, R'') 
minimizes the functional 



J"l{u,v,R) = I^^'^^'^iu.v) - R33 / ( / f^dxX{x')dx' +r{K), 

(3.89) 

where is defined in l3Jl\) . 

3. The minimum of the functional J*^ i.e. exists in the space VF^'^(w;M^) x 

W^'^^^oj) X S0(3). If e W^'^{n'^]'R^) is a minimizing sequence (not 
1 



relabeled) of -^J^ then we have that 



\UB. -^J^{y^) = min J^{u, v,Ii), if h-'^E^ ^ I, 

h^O E'^ uGW1.2(aj;R2), veW^-^{uj), ReSO(3) 

(3.90) 

lim ^ J''(y^) = min j'^ (u,v,R), if h-'^E^' ^ 0. 

h^O E'^ iiGW'i.2(w;M2)_ ,)eW2.2(a;), ReSO(3) 

(3.91) 

Proof. The proof goes in the same direction as the proof of Theorem 2.5. in 
[15] . If we take y'' = 'RP'x where R'* is chosen such that 



h Ich 



f^,{x){R''xhdx, 



we see the left inequahty in (I3.86p . Let us now take the minimizing sequence 
y^. Using Theorem ([3|) and the coercivity property of energy density function 
we find R'' e SO (3) such that 



KVy") R^|U.(,,^„3x3) < Ch-'JV^{yh). (3.92) 
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Taking — f^^ y'^ o 0^ o P'* we conclude, by using Poincare inequality the 
boundedness of V©'', {V0'')-\ VP'* that 



(R'')V-c''- id\\ 



L2(r>:R3) 



Using p.80p we conclude 



< C^m'^fiy'^ o 0'^ o P>^) -c'^-0^0 P"|U.(^2;K3) 

< C^m'^fViy'^ O 0'^ O P") - V(0'^ O P")|U.(^2;K3X3) 

< C-^||(R'')^(Vy'') o 0" o P'^ - I|U.(r,;M3x3) 
<C^h-\flNy^). (3.93) 



^ ^ (R")^ I I .xdx- ^. (3.94) 



h 



Using p.80p , p.87p , the left inequality in p. 861) and the definition of m'* we 
conclude that there exists C > such that 



l'\y'') < C{E'' + ^/e^^JJh^)). (3.95) 

From this we conclude that there exists C > such that I^{y^) < CE'^. 
Using (|3.93p we conclude -^J^iy^) > — C. From p. 871) we conclude the right 
inequality in (I3.86p . 

Everything else is the consequence of the fact that -^J^ ^ J° (see [TUl 

E]). 

To prove the lower bound we take y^ e H/^'^ (!?'»; R^) such that 
liminf/i^o ■^J'^iu^) < +00. We take subsequence (not relabeled) such that 
liuih^o j'^ {y'^) = liminf/i^o -^•^'^ (?/'') • In the same way as in p.94p we 
conclude I'^iy'^) < CE'^. 

Using Lemma m we find R'' e S0(3), c'' e M, m e W^''^{Q-M^), v G 
W'^''^{u)) such that for the subsequence (not relabeled) y^ := {ti^)'^y^ — c'', 
u^, v'^ it is valid u'^ it in VF^'^(a;;R^) and ^ v in W^'^(a;). Let us take 
subsequence (also not relabeled) such that R'' R. Let us write 
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fUx){Ii''xhdx 



(3.96) 
(3.97) 



Let us prove that 
1 



1/2 



1/2 



f3{x',X3)dx3)v{x')dx'. (3.98) 



Let us consider 
1 



f3ix)iy3i^) ~ x)dx 



1/2 



-1/2 



f3{x' , X3)dx3)v{x')dx' 



< 



{f^o0^^opf^)ix) h 



hVE^ 



{flo0^op^){x: 



((y^0"oP'')(:r)-(0^P")(x)) 
-v(x')J det ({V0''){P^{x))^ dx 
v{x') det ({V&^){P''{x))] dx - 



hVE'^ 

f3{x)v{x') det [{V0^){P\x))) I + C{f^) 



(3.99) 



The first term in p. 991) goes to by Lemma [31 the boundedness of 

and the wealc convergence in p. 801) and the second term goes to by the 

weak convergence in (|3.80l) . Thus we have proved p. 981) . From the fact that 

\^j'^{y^)\ is bounded we conclude that -f^{m^^ — f^-^^ f'^{x)(R'^x)3dx^ is 

bounded. From that we conchidc R e A^. Using Corollary Q we conclude 

itJ?f P'^y"^ ^ itJ?f P'^y'^ + iTo 7^ir miy'.i^) - -3)rf- 



:tdx 



> 



J°(m,w,R), if /i-^e:'' ^ 1 
J^(m,w,R), if /i-'^i;'' ^ 



To prove the upper bound for u e W^'^{uj;M.'^), v G VK^'^(cj;] 
choose y'^ from Theorem 2] and R'' — > R such that 



R e we 



r(R) = lim 



1 



h^o hE^ 



/5(x)(R''x)3dx 



(3.100) 



We define j/'' = R''y'\ We decompose the total energy functional 




(y (a;) — x)dx 

(3.101) 
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It is easily seen, by using Lemma 121 that 




From the decomposition p.lOll) we conclude, by using (jS.lOOp and (j3.103p . 



that ^J'^iy'') J"(m,w,R), if h-^E'' 1 and -^J^{y'^) JO(m,?;,R), 



if h~^E^ — 0. Thus we have proved the P-convergence result. Existence of 
the minimizer follows from the r'-convergence theory, but can also be proven 
independently (see Lemma[S]and the references in [I] Section 5.12]). 

Remark 7 The analysis in the situation when we have also in-plane forces 
(dead loads) is more complicated. The problems that appear in these situations 
are similar to the ones observed for the classical plate in [19]. Namely, if we 
want to add in-plane forces, they should be of order in all the situations 
h~^E^ — >• oo or h~'^E^ — >■ 1, since u is of order j^. From the expression p.94p 
we can only conclude that I^{y^) < Ch^. Thus we do not know which model 
is appropriate for these forces and we would have to impose some stability 
condition (see [T^). The situation is better in the case h~'^E'^ — >■ (the order 
of the in-plane forces should then be V E^), despite the fact that we can not 
conclude directly the right order of the internal energy from the expression 
p.94p . This is because the body can perform large rotation and mix the in- 
plane and the normal forces. This can be prevented by imposing appropriate 
boundary conditions (see |19)1. 

Lemma 5 For /g e L'^{Q) such that J^^ f^^O and g : S0(3) MlJI+oo} 
bounded from below, lower semi- continuous and not identically equal to {+00} 
on the set Ai'^, the Junctionals 



have the minimizers in the space W^'^{u]\^^) x W'^'^{u]) x AIJ,, where Ai'^ is 
any closed subset of the set 




(3.105) 




VF e so(3)}. (3.106) 
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Proof. Sequentially weakly lower semi-continuity of the functionals is the 
direct consequence of the convexity of the form Q2 and Rellich-Kondrachov 
compactness embedding theorems (see e.g. |22j). The boundedness of the min- 
imizing sequence is the only fact that has to be proved. We shall do it for 
the functional J". Let us take the minimizing sequence (it^, R'*) such that 
J"(u.'', w'*, R'*) — > inf < -foo. Since Q2 is positive definite on symmetric 
matrices, we have that I^'^""^ {u,v) > Ci||(V')^t;||L2(f2). By the Poincare in- 
equality there exist a'^ G R^,b'' e R such that - a'' • x' - b''||vF2,2(f2) < 
C'2||(V')^w''|ji2(j7). Let us note that the condition in the definition of the space 
p. 1061) for Q33 7^ implies f^^i — and f^^i — 0. Thus in either case 
we have R33 f^{x'){a''^ ■ x' + h''')dx' — 0. From this we can easily deduce 
(by the standard argumentation) that \\v^' — ■ x' — b''|| 14/2,2(17) < C. The 
problem is, however, that the functions of the type a'' • x' affect the part of the 
functional (52(i(V'M'' -I- (V'tt'')^ + V'v^ (g) V'u'' + Vu'' (g) V'O + V'O ® V'v^)). 
But we have 

v'm'* + {Vu'^f + Vv^ (g) Vv^ + Vv^ ® ve + V9 ® v'w'' = 
v'lt^ + v'(m^)^ + Vv^^ ® Vv^^ + v'w^ (g) ve + \/'e ® v'w^ 

where 

u'l = u'' + {v'' + e)a'' + ^('^'^ • = v'' -a^ -x - b'\ 

Thus we conclude that 

J"(tt^, R'') = J"(M^ R'') ^ inf J°, 

for ||w^||vK2,2(r2) < C. Now we take Aj? e M2x2^ y^h ^ ^ ^^^^ ^^^^ ^^^^ 
Korn's inequality 

- A,^x' - bJ|k-(r2;M^) < C3II sym(V'«,^)|U2(^^R2). (3.107) 

If we define w^'^ = u'^ — A^x' — b^ we have 

J°(m;?„w^R'') JO(M'',^;^R'') ^ inf J°, 

for ||w''||vi/2,2(j-2) < C. From the boundedness of and the fact that 

/^^"^(m^„«^?,R'') > Ci\\Vu'^, + V{u'^,f + Vv'^®Vvl' + 

+v'v'^ ® ve + ve ® v'«,^||L2(fi;R2), 

we easily conclude the boundedness of IIV'mJ?^ + V'(ttJ?j,)'^||i2(f^2;R2) and from 
p. 1071) the boundedness of |jM^^||iyi,2(f2;a2). From this it follows that the se- 
quence {u'^^,v^ ,'Rp-) has its weak limit. Due to the sequentially weakly lower 
semi-continuity of the functional J", the weak limit of the minimizing sequence 
is a minimizer. 
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Remark 8 If we want to include tangential forces and displacement boundary 
conditions, then we could obtain existence result for the functional J° only for 
small tangential forces (see [8ll9ll25)). This is not the case for the functional j£ 
where we can, by using the generalized Korn's inequality (see [T]), conclude the 
coercivity of the functional in the appropriate space including the displacement 
boundary conditions. 
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